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Abstract. We determine the behavior of the general solution, small or large, 
of nonlinear first order ODEs in a neighborhood of an irregular singular point 
chosen to be infinity. We show that the solutions can be controlled in a ramified 
neighborhood of infinity using a finite set of asymptotic constants of motion; 
the asymptotic formulas can be calculated to any order by quadratures. These 
constants of motion enable us to obtain qualitative and accurate quantitative 
information on the solutions in a neighborhood of infinity, as well as to deter- 
mine the position of their singularities. We discuss how the method extends 
to higher order equations. There are some conceptual similarities with a KAM 
approach, and we discuss this briefly. 

1. Introduction 

The point at infinity is most often an irregular singular point for equations arising 
in applications Q Within this class of equations, there are essentially two types for 
which a global description of solutions exists: linear systems and integrable ones. 
However, in a stricter sense, even for some linear problems global questions such 
as explicit values of connection coefficients are still open. The behavior of the 
general solutions of linear ODEs has been thoroughly analyzed starting in the late 
19th century, see [24] and [38] and references therein. After the pioneering work of 
Ecalle, Ramis, Sibuya and others the description of their solutions in C is by now 
quite well understood [23l [22l [H El HUl EH Ell US] . 

Integrable systems provide another important class of systems allowing for global 
description of solutions. The ensemble of integrable systems is a zero measure set 
in the parameter space of general equations: a generic small perturbation of an 
integrable system destroys integrability. Nonetheless, integrable equations occur 
remarkably often in many areas of mathematics, such as orthogonal polynomials, 
the analysis of the Riemann-zeta function, random matrix theory, self-similar so- 
lutions of integrable PDEs and combinatorics, cf. [8], [18] -[20], p] [TTJ EE21 O [20] , 
@D] • However, even in integrable systems, achieving global control of solutions 
in a practical way is a challenging task, and it is one of the important aims of the 
emerging Painleve project [9]. 

In nonintegrable systems, particularly near irregular singularities, our under- 
standing is much more limited. Small solutions are given by generalized Borel sum- 
mable transseries; this was discovered by Ecalle in the 1980s and proved rigorously 



A singular point of an equation is irregular if, for small solutions, the linearization is not of 
Frobenius type. By a small solution we mean one that tends to zero in some direction after simple 
changes of coordinates. 
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in many contexts subsequently. Transseries are essentially formal multiseries in 
powers of l/x k ' e" A ^, and possibly x' 1 \ogx; see again [2311221301231231311112] 
and [T3] • Here x is the independent variable and are eigenvalues of the lin- 
earization with the property Re (Xjx) > 0. In general, only small solutions are well 
understood. However, for generic nonlinear systems of higher order, small solutions 
form lower dimensional manifolds in the space of all solutions, see, e.g., [161 . The 
present understanding of general nonlinear equations is thus quite limited. 

We introduce a new line of approach, combining ideas from generalized Borel 
summability and KAM theory (see, e.g. [3) for the analysis near infinity, chosen to 
be an irregular singular point, of solutions of relatively general differential equations 
with meromorphic coefficients. Applying the method does not require knowledge 
of Borel summability, transseries or KAM theory. 

For small solutions, in |17) it was shown that in a region adjacent to the sector 
where the solution, y, is small, y{x) is almost periodic. In this sense y becomes 
an approximately cyclic variable. In the x-complex plane, the singular points of 
y are arranged in quasi-periodic arrays as well. The analysis in |17j covers an 
angularly small region beyond the sector where y is small. Looking directly at 
the asymptotics of y beyond this region would require a multiscale approach: y 
has a periodic behavior-the fast scale, with 0(1/ x) changes in the quasi-period. 
Multiscale analysis is usually a quite involved procedure (see, e.g., [7]). 

It is natural to make a hodograph transformation in which the dependent and in- 
dependent variables are switched. As mentioned above, in the "nontrivial" regions, 
the dependent variable is an almost cyclic one. The setting becomes somewhat 
similar to a KAM one: there is an underlying completely integrable system, and 
one looks for persistence of invariant tori. Adiabatic invariants are simply the con- 
served quantities associated with these tori. Evidently there are many differences 
between the ODE setting and the KAM one, for instance the fact that the small 
parameter is "internal", l/x. 

In this work we restrict the analysis to first order equations, mainly to ensure 
a transparent and concrete analysis. In theory however, the method generalizes to 
equations of any order, and we touch on these issues at the end of the paper. 

We look at equations which, after normalization, are of the form dy/dx = F(z, y), 
z = l/x. with g bi-analytic at and F y (0) = 1. 

We show that in any sector on Riemann surfaces towards infinity, the general so- 
lution is represented by transseries and/ or, in an implicit form, by some constant of 
motion. In fact, on large circles around x = 0, the solution cycles among transseries 
representations and ones in which constants of motion describe it accurately. The 
regions where these behaviors occur overlap slightly to allow for asymptotic match- 
ing (cf. Corollary [5]) ■ The connection between the large x behavior and the initial 
condition is relatively easy to obtain. 

Let /? = F zy (0). The constants of motion have asymptotic expansions of the 
form 

(1) C(x,y)~x-/3logx + F {y)+x~ 1 F 1 (y) + --- + x- : >F j (y) + --- , x -> oo. 

Clearly, under the assumptions above, the solution y can be obtained asymptotically 
from ([I} and the implicit function theorem. The requirement that C is to leading 
order of the form fi(x) + f2{y), determines C up to trivial transformations, see 
Theorem [5] and Note El 
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The functions Fj are shown in the proof of Theorem [5] to solve first order au- 
tonomous ODEs, and thus they can always be calculated by quadratures. 
To illustrate this, we use a nonintegrable Abel equation, 

(2) v! = — t. 

We note that there is no consensus on how nonintegrabilty should be defined; 
for |2]), it is the case that the equation passes no criterion of integrability, including 
the poly-Painleve test, and that there are no solutions known, explicit or coming 
from, say, some associated Riemann-Hilbert reformulation. 

The Abel equation has the normal form (see 21 where further details about this 
example are given) 

(3) y,+Sy3 -l + h y = °- 

Regions of smallness are those for which y approaches a root of 3j/ 3 — 1/9; in these 
regions, y is given by a transseries [17] . Otherwise, y has an implicit representation 
of the form 

(4) y = i exp ^ - C - x + i log x + - -jr^J arctan ~ ^ 

- - » + 5 ^( 9 » 2 + % + '» + ; (^Wj + k + 0(im) ) + i, 

obtained by inverting an appropriate constant of motion C (see (|39p ); for the values 
c£0,F o ,Fi see S3 

While in a numerical approach to calculating solutions the precision deteriorates 
as x becomes large, the accuracy of ([T]) instead, increases. In examples, even when 
(JXJ) is truncated to two orders, ([T]) is strikingly close to the actual solution even for 
relatively small values of the independent variable, see e.g. Figure 0J 

The procedure allows for a convenient way to link initial conditions to global 
asymptotic behavior, see e.g. (|4"Tj) . 



1.1. Solvability versus integrability. First order equations for which the asso- 
ciated second order autonomous system is Hamiltonian are in particular integrable. 
Indeed, by their definition, there is a globally defined smooth H with the property 
that i^f- + y^- = 0, that is H(x(t), y(t)) = const, providing a closed form implicit, 
global representation of y. While the differential equation provides "infinitesimal" 
information, H -effectively an integral- provides a global one. 

Conversely, clearly, if there exists an implicit solution of the equation or indeed a 
smooth enough conserved quantity, the equation comes from a Hamiltonian system. 

What we provide is a finite set of matching conserved quantities, analogous to 
an atlas of overlapping maps projecting the differential field onto the trivial one, 
H' = 0. They give, in a sense, a foliation of the phase space allowing for global 
control of solutions. With obvious adaptations, this picture extends to higher order 
systems. In integrable systems there is just one single-valued map and the field is 
globally rectifiable. In general, the conserved quantities may be branched and not 
globally defined. 
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1.2. Normalization and definitions. Many equations of the form y' = F(y, 1/x) 
with F analytic for small y and small 1/x can be brought to the normal form 
y' = Po(y) + Q(y, 1/x) by systematic changes of variables, see e.g. [16] . |15j . 

The assumptions are that Q(y, z) is entire in y and analytic in z for small z, and 
0(y 2 ,yz 2 , z) for small y and z and that Pq is a polynomial. We assume that the 
roots of Pq are simple. It will be seen from the analysis that a more general Pq can 
be accommodated. We thus write the equation as 

(5) y' = £ ^ = Qx(y, l/x) = P (y) + Q(y, 1/x). 

k=0 

Definition 1. • A formal constant of motion of |5]) for x — > oo in an unbounded 
domain PcC 2 or on a Riemann surface covering it, and in which to leading order 
in 1 fx the variables x and y are separated additively is a formal series 

(6) C(y,x)=A(x)+F (y) ' Fl{v) < > Fj{y) < --- 
such that we have 



±C(y(x),x) = 0(x-°°) 
in the sense that, for any j , Fj and Hj defined by 

(7) ^±#^ := A\x) + D x ( F (y) + F ^ y) ■ ■ ^ 



x 1+L \ x x 3 

are uniformly bounded in T>; here D x is the derivative along the field, 
D x F(x, y) = VF- (1, Qi) = F x {x, y) + F y (x, y)Q 1 (y, 1/x). 
See also (|10[) below. 

• An actual constant of motion associated to C in D C C 2 is a function C so 
that C(y,x) ~ C(y,x) as x — > oo and -4-C(y(x),x) = for all solutions in T>. 

Note 2. It will be seen that there is rigidity in the form of the constant of motion: 
if the variables in C are, to leading order, separated additively as in ([6]), then, up 
to trivial transformations, we must have 

(8) A(x) = — x + a log x 

where a is the same as the one in the transseries expansion of the solution, see 
Proposition 



1.3. Finding the terms in the expansion of C. Using ((HJ and truncating (O 
at an arbitrary n > 2, let 



Fk{y) 

x k 



(9) C n {y,x) =: _ x + a log x + F (y) + ^ ^ 

fc=i 

We can check that D x C n satisfies 

(10) D x C n = -l + P F> + a + PlF ° + P ° Fi 

x 

" (1 - k)F k ^ + YlU Pk-jFj -nF n + J2U EZo Pn+k+i-jF^x- 

k=2 
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(cf. (|5|l) where the numerator of the last term is H n by definition. In order for C to 
be a formal constant of motion, the coefficients of x~ 3 , j = 0, 1, 2, . . . must vanish, 
giving 

(11) Fo(v) ' 



(12) F[{y) = 



Po(y) 

a + Pi(y)F£(y) 



My) 



It follows in particular that Fq ^ and i*o is bounded in T>. In solving the differen- 
tial system, the constants of integration are chosen so that are indeed uniformly 
bounded in y, see (|25l) . 



1.4. Solving for y(a;). The expression C n is an approximate constant of motion ; 
we thus can find an approximate solution y n by fixing C n — K. We then write 

(14) G(y; K) := F (y) - K - x + a \ogx + ^ = 

fc=i 

and we note that in the domain relevant to us (Si , see Theorem [4] below) the 
analytic implicit function theorem applies since 

dG 1 1 , . 
15 — = —— + -E 1 (y,x) 

oy Po(y) x 

where Pq is away from in our domain, and for some E\ which is bounded in T> 
by ((9J since Y is bounded. Writing y = y n in ([M)) and using the analytic implicit 
function theorem, treating 1/x as a small parameter, we get 

n*\ rt ^ , G x (x;K) , g^^jQ , ffn(x;iC) 
(16) y n = G (x;K) + + ••• + — + ^ 

where and the G^'s are bounded. In the same way it is checked that y n 
is solution of ([5]) up to corrections R n (x; K )/x n+1 , that is, y' n — Qi(y n ,l/x) = 
—R n (x;K)/x n+1 where R n is bounded. 
Let pi, . . . ,p m be the distinct roots of P . 

Let TZ y be the universal cover of Y = C\{pi, ...,p m }. Let it : X — > Y be the 
covering map. 

Definition 3. • An elementary y-path of type 

a = (ai, ...,a m ,a m+ i, . . . ,a mfc ) € Z mk ,k e N 

is a piecewise smooth curve 7 in TZ y whose image under ir turns a\ times around 
pi, then «2 times around P2, and so on, a m times around p m , then again a m+ i 
times around p\ , etc. Note that a is in fact an element of the fundamental group. 

• A y-path of type a is a smooth curve 7 obtained as an arbitrary forward 
concatenation of elementary y-paths of type a. More precisely, a y-path of type 
a is a map 7 : [0,oo) — > 1Z V so that, for any N e Z + , 7|rjv,jv+i] i s an elementary 
y-path of type a. We will naturally denote by 7|[o, a ] subarcs of 7. We see that 
y-paths are compositions of closed loops in the complex y domain. 
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• S r is a regular domain of type a or an i?-domain of type a, if it is an 
unbounded open subset of 1Z V that contains only images of y-paths of type a. Thus 
the image of any unbounded y-path of type a' 7^ a is not a subset of S r . 

Note. In our results we only need y-paths with the additional property that 
x {y) ~* 00 along the path. 

To take a trivial illustration, in the equation y' = y an example of a y-path along 
which x — > 00 is t n- exp(it), t > 0. 



2. Main results 

2.1. Existence of formal constants of motion. Under the assumptions at the 
beginning of ^1.21 we have 

Theorem 4. Let S y be an R-domain of type a, and 

Si = {y e S y : \ir(y)\ < M Q and \n(y) - p k \ > £ for all k} 

where Mq > is an arbitrary constant. Let C be the union of m circular paths 
surrounding a^eZ times the root pk, k = 1, . . . , m, chosen so that 

(!7) / -BT^dV ± 0. 



JC My) 

Then, if Rq is large enough, there exists a formal constant of motion in 

V 1 = {(x,y):\x\>R ,y€S 1 } 

of the form ([6]) . The terms F k in the expansion of C in ^ can be calculated by 
quadratures. 

Actual constants of motion are obtained in Theorem [5j 

Consider now a set S of curves 7, |7(i)| — > 00 as t — > 00, with the property that 
for all ti < ti and all n (which is in fact equivalent to for n — 0, 1 ) 



(18) 



7(*2) Q 

Re / — Qi(y,7(*))l 1/ =3/„(7(t))7 / (*)* 

7(*i) °y 



^61og(|7(f 2 )/ 7 (ii)|+l), 



where b > is a constant, and such that there is an M so that for all n we have 
\y n \ < M along 7. Here M can be chosen large if x is large. Note that S contains 
the curves 7(f) so that y n (j{t)) is an a-path. Indeed, by p^|) , in this case, the 
integrand in (fT8| is of the form p°[^j dy + (9(1) and hence the integral equals 
2niN + 0(log(| N\ + 1)) for large N where N is the number of loops. 

Theorem 5. Assume C in ^ is a formal constant of motion in a region T> — 
SnDi. Then there exists an actual constant of motion C = C(x,y) defined in the 
same region, so that C ~ C as x — > 00. 

2.2. Regions where Pq(u) is small. Assume x is large and \P (y(xo))\ < e is 
sufficiently small. This means that for some root ru of Pq we have \y(xo) — r^\ < ei 
where t\ is also small. Without loss of generality we can assume that r k = and 
x £ R + since the change of variables yi = y — r k , x — xie 1 ^ does not change the 
form of the equation. Assume also that after normalization the stability condition 
ReP^O) < holds. Again without loss of generality, by taking y^ — ayi we 
can arrange that Pq(0) = — 1. The new function Q in ([5]) will have the form 
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y 2 Qi(y, l/x) + x 2 Q2(y, where Qi and Q 2 are analytic for small y and 1/z. 
As a result, the normalized equation assumes the form 

(19) y' = -y + fo(x) + ^ + y 2 Q 1 {y, l/x) + x~ 2 Q 2 (y, 1/x). 

x 

We also arrange that /o = 0(x~ M ) as x — > 00, for suitably large M; this is possible 
through a change of variables of the form y 2 = yz + J2k=i c k x ~ k , where the c k s 
are the coefficients of the formal power series solution for small y. 

Proposition 6. [see [TB] Theorem 3] Any solution of (flT)]) that is o(l) as x — >• 00 
along some ray in the right half plane can be written as a Borel summed transseries, 
that is 

00 

(20) y(x) =Y / C k x ka+1 e~ kx y k 

k=Q 

where y k are generalized Borel sums of their asymptotic series, and the decomposi- 
tion is unique. There exist bounds, uniform in n and x, of the form \y n {x)\ < A k , 
and thereby the sum converges uniformly in a region 1Z that contains any sector 
S c := {x : I arg x\ < c < ir/2}. Note that Theorem 3 in [16] applies to general n-th 
order ODEs. 

Proposition 7. (i) If, after the normalization above, y(xo) is small (estimates can 
be obtained from the proof), then y is given by (|20[) . 

(ii) C{y{x),x), obtained by inversion of &20\) for large x in the right half plane 
and small y, is a constant of motion defined for all solutions for which y(xo) * s 
small (cf. (i)). 

Proof, (i) We write the differential equation in the equivalent integral form 

(21) y = F Q (x)+y e~( x -*°\x/xo) a 

+ e~ x x a [ e s s' a [y 2 (s)Qi(y(s), l/s) + S - 2 Q 2 (y(s), 1/s)] ds, 

Jx 

where Fq(x) = 0(x~ M ) (M can be chosen arbitrarily large in the normalization 
process, [TB]) and Fo(a;o) = 0. It is straightforward to show that for (|2"TT) is con- 
tractive in the norm \\y\\ = sup xeSc \x M ~ 1 y(x)\ (see the beginning of this section) 
and thus it has a unique solution in this space. Hence, by uniqueness, the solution 
of the ODE with y(xo) = yo, has the property y(x) — > as x — > 00. The rest of (i) 
now follows from |16j . 

(ii) We see from Proposition [5] that y(x; C) is analytic in a domain of the form 
S c x O p (As usual, D p denotes the disk of radius p.) We look at the rhs of (|2"0"|) as 
a function H(x, C). It follows from [16] that y\{x) — a; _1 (l + o(l/x)). By uniform 
convergence, we clearly have 

( 22 ) -qc = ^ kC k - x x ka + x e- kx y k = e~ x x a {\ + o(l)) + 0. 

k=0 

The rest follows from the implicit function theorem. □ 
As a result of Theorem [S] and Proposition [7] we have the following: 

Corollary 8. // Go in (|16|) approaches a root of Pq and x is large enough, then y 
enters a transseries region, where the new constant is given, after normalization, 
by Proposition^ (ii); thus the constants of motion in different regions match. 
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3. Proofs and further results 

3.1. Proof of Theorem^ Let (xq, yo) € T>\. Recalling (fTO]) . we see that (fTTj) has 
the solution 

rv i 

Mv) = / 77TT ds + c o 
Jyo P o(s) 

(we take Co = since it can be absorbed into the constant of motion). Eq. (|T2|) 
gives 

(23) F X (y) = /^y) + d := - / - f'\ (s) ds + Cl) 



v a + 



where to ensure boundedness of Fi(y) as the number of loops — » oo, we let 



and ci is determined to ensure boundedness of F2 (cf. (|25p). Inductively we have 
(24) 



Fk+i(y) = / =-7-^ + c fc+1 =: fk+i{y) + Cfc+i 



for 2 < k + 1 < n, and, to ensure boundedness of Fk+i(y) as the number of loops 
— > 00 we need to choose 

. -kfk+T,U p w-j(y)fj(y)^ 

(25) c k = W 

kJc m dy 

for 1 < k < n — 1. 

It is clear by induction that every singularity of F k (y) is a root of Pq . To complete 
the proof we need to show that the Fk 's are bounded in T)\ : 

Lemma 9. Assume y G Si. For deg(Po) > 1 owd 1 < k < n we have 

\F^(y)\<k\ and \F k (y)\ < k\(\y\ + 1) 

where, as usual, < means < up to an irrelevant multiplicative constant. 

Proof. We prove the lemma by induction on fc. Note that in (|2"3"|) and (|2"4"|) the 
integration paths can be decomposed into finitely many circular loops C and a ray, 
slightly deformed around possible singularities, which implies 

\Fi(y)\ <log|y| + l<|y| + l 

and 

l^-(y)l< [ V \FL(s)\ds 

Jyo 

where the integration path is a straight line (possibly bent as above) 
We see from (EH) that 



< (fc-i)|F fc -i(y)l y }| < (fc-i)!^-!^)! t4 
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The conclusion then follows by induction. Note that the the last term of (fT0|) 
satisfies 



-nF n +J2J2 Pn 

+k+l-jF' J X 

j=0 k=0 



< 



(n + l)!(|y| + l)|P (y)| 



□ 



3.2. Proof of Theorem [5]. Let y(x;K) = y n (x;K) + 5(x;K), where y n is given 
in (|TB)) . We seek S so that y is an exact solution of ([3]) in T>. 

Let (j>(y, 8) be the polynomial satisfying Qi(y + S,x) — Qi{y, x) = Qi^ y (y, x)S + 
S 2 (j)(y,S,x) where Qi, y (y,x) := 9Ql ^' x) . We obtain 

(26) 6 6= — — \-(p(y n ,5,x)5 =: E(x; S(x); K), 

x ay x n+i x 

where R =: R n is defined after ([TBI) ; both R and <f> are, by assumption, bounded. 
In integral form, (|26|) reads 

(27) 5(x) = f X ^ef* Q^Vn(t),t)dt E ^ s . g^. K ^ ds 



where the integrals are taken along curves in D. Using (fT5|) we see that (|2~T1) is 
contractive in the norm ||<J|| = sup |x|™|(5(x)| in an arbitrarily large ball, if 

|x|>|ai|;x£Z) 

|a;i| is large enough and n > b2 b+1 . 

Thus ([2T]) has a unique solution and, of course, S(x) is the limit of the Picard 
like iteration 

(28 ) 6 = r ^ef*^y^ dt ^£±ds 



5i= I ^-efs Q ^^W dt E(s;6 (s);K)ds 



J' 

etc. 

By (|2T|) (5 is a smooth function depending on (x, K) only, and 8 = 0(x~ n ). Smooth- 
ness is shown as usual by bootstrapping the integral representation (j27l) . 

Now we have, by (fl5)) . dxynix; K) — P (y n )(l + o(l)). We can easily check that 
8k8(x, K) = 0(x~ n ). This is done using essentially the same arguments employed 
to check contractivity of the integral equation for 5 in the equation in variations for 
Sk, derived by differentiating (f2l)|) with respect to K. We use the implicit function 
theorem to solve for K, giving K = K(x,y), a smooth function of (x,y). It has 
the following properties: K(x,y(xj) is by construction constant along admissible 
trajectories and by straightforward verification, i.e. comparing K with C, we see 
that it is asymptotic to C up to 0(x~ n ). It is known that if a function differs from 
the nth truncate of its series by 0(x~ n ) for large n, then in fact the difference is 
o{x~ n ) (cf. [U] Proposition 1.13 (iii)). 

3.3. Position of singularities of the solution. It is convenient to introduce 
constants of motion specific to singular regions; they provide a practical way to 
determine the position of singularities, to all orders. 



10 



O. COSTIN, M. HUANG AND F. FAUVET 



Definition 10. We define a simple singular solution path 7(5) : [0, 1) — » 1Z y 

to be a piecewise smooth curve whose projection 7r(7[0, 1)) € C is unbounded but 
turns around every pk only finitely many times. 

A simple singular solution domain S s is the homotopy class of any simple 
singular solution path, in the sense that any two unbounded paths in S s can be 
continuously deformed into each other without passing through any pk ■ 

Proposition 11. Let mo = deg(Po) > 2, S s be a simple singular solution domain, 
and T>2 = {(x, y) : \x\ > R,y G S s , and \y — Pk\ > e for all k}. Assume that 

\3M < u-« 

\Po(y)\~ m 

for large y, for some q > and all k > 1. Note that this needs only be true in S s , 
which could be an angular region. 

Then there exists in T>i a formal constant of motion of the form 

(29) c(y,x)=x + F (y) + ^ + --- + ^ + --- , 

x x J 

where Fk(y) are single valued as y — > 00. Furthermore, any simple singular solution 
path passing through some arbitrary {xo,yo) tends to a singularity, whose position 

X sitiq S(XtlsJ~L€S 

(30) x sing = C n (y Q ,x ) + O 



sing — ^nVi/U? ^U; 1 \ n+1 

for all neFt, where C n is C truncated to x~ n . 

Moreover, if there are only finitely many nonzero Pk, then there exists in T>i a 
true constant of motion of the form \29\) . i.e. the sum is convergent for large \x\. 

Proof. The proof is similar to that of Theorem SJ 

In order for C to be a formal constant of motion, we must have 

(3i) H(v) = - ' 



Poiv) 



(32) m _ *-^V-T£^WM {lstin , 

Po(y) 

We solve successively for the F^ and obtain 

rv 1 

My) = / -ttt-Js 

J 00 Pa{s) 

where the integration path lies in S s . Clearly Fq is bounded and single valued as 
y — > 00. 

Inductively we have 



(k-l)F k _ 1 ( 8 )-Y^Pk- 3 <8)F>J8) 



(33) F k (y)= J - ' ds 

for 1 < k < n. 

To prove the rest of the proposition, we need the following lemma: 
Lemma 12. Assume that y G S s . For 1 < k < n we have 

k\ 



\n(v)\ 



< 



V 



mo+q 
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\F k (y)\ 



< 



m +q-l 



\v 



as y — > oo . 

Furthermore, if Pk = for k > ko > 0, then 

\K(y)\ < — 



w 



mo+q 



\My)\ 



< 



|y|mo+g , -l ' 

Proof. The estimates are obtained by induction on k. Note that (|32|) implies 

fc-i 

I I ■ . 111/', lyil 

l^(y)l 



.(MM 



I. v I 



provided that the assumptions of the lemma hold for 1 < j < k — 1. 

If Pfe = for k > ko > 0, we again show the lemma by induction. 
Assume that for < / < k we have 



2+1 



\Fl(y)\<(colo) l J2{-_ 1 



(this is obviously true for I = 1). 
This implies 



i=i 



fe+i 



-l+.?( 1 - m o)-<? 



|^.( 2 /)|<( Co A ; o) fc E( ? „ 1 
Thus it follows from ([32]) that 



\j(l-m )-q 



^ vi - 1/ j( m o - 1) + 9 



fc^fc(y) - S?= ma x{fc-fco+i,o} p k+i-j(y)Fj{y) 
My) 



where, by the induction assumption, the first term satisfies the estimate 



(34) 



kF k (y) 



Po(y) 



< c\k 



F k (y) 



k+2 



< r k k k+1 n V — - 



(fc + D(A) 



mo)— 9 



(j-l)(mo-l)' 



fc+2 /, , 

<« +1 C iE(- 

i=2 v 

where cq > 1 + ci . Note that the last inequality follows from 



-l+j(l-m )-g 



(fc + 1) 



fc + 1 

i-i 



0' - !)• 



A- 

VJ - 2y 

The second term is easy to estimate, since it is clearly bounded by 

fe+i 

-m a )-q 
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Since c\ is fixed, we can assume that cq > 1 + c\ , and we have 

fc+2 



i^ +1 (2/)i<( C ofc ) fc+i E( fc !!)i^ 1+ ^ 

j=i V / 



m )—q 



This shows the second part of the lemma. □ 
Now since 

-nF n + Y,UY;Zo p n+k + i-jFlx- k < |P ( y )| 



\D x C n 



r n+l 



|a:|n+l|y|mD+« 



(cf. UOp , the estimate for x S i ng follows immediately from integrating D x C n from xq 
to x a i n g — C n (oo,x s i ng ) along the simple singular solution path. 

□ 



<\y\~ q 



Remark 13. The condition 

\P k (y)\ 
\Po(y)\ 

is not the most general one for which there exists a formal constant of motion in 
a simple singular domain. However, this condition is frequently satisfied by ODEs 
that occur in applications (see <Q). In such cases we can easily use \30\) to find the 
position of the singularity (see e.g. \4-l\j ). 



4. Example: the nonintegrable Abel equation J2]) 

To illustrate how to obtain information of the solution of a first order ODE 
using Theorem[3]and PropositionQTJ we take as an example the nonintegrable Abel 
equation ©. Normalization is achieved by the transformation x = — (9/5)A 2 t 5 ^ 3 , 
A 3 = 1, u{t) = A 3 l b (-Uh) 1 l b x 1 l b y(x) PI], yielding 

(35) y> = -Zy* + \-y 

9 5a; 

Obviously (f3"5|) satisfies the assumptions in Theorem 0] and [Til with P (y) = 

-3y 3 + i and Pi(y) = -|. 
9 5 

1 (_l)2/3 (-1) 4 / 3 ,_, 

The three roots of Po are -, , and . It is known [17] that there 

exists a solution in the right half plane H that goes to the root — as 1 — > 00. 

Similarly, there are solutions that go to the other two roots in other regions, which 
we will explore in £14.31 In those cases, the behavior of the solution follows from 
Proposition [6] (see also [T7]). However, there are also solutions that do not go to 
any of the three roots. In these cases, the formal constant of motion will be a useful 
tool to describe quantitatively the behavior of the solution. 

4.1. Constants of motion in i?-domains. (cf. Definition [3]) . First we choose an 
elementary solution path along which the solution y to (135|) turns around the root 

- clockwise as shown in Fig [JJ and [2J 
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FIGURE 1. Solution y{x) with y 
from l+5i to 1.5+5CK to 1.6+120i 



= 1.1 along the line segments 




Figure 2. Real and imaginary parts of y(x). The upper curve 
is the real part, the lower curve is the imaginary part, and the 
straight line is the root 1/3. 

For simplicity we calculate the first two terms of the expansion ([9]). We have 
(36) 



F Q (y) 
(37) 

a 

(38) 
Fi(y) 



- jdy = \/3arctan 

-3lf + n 



fc 5 (-3y3+±)2 d y _ 1 
~ 5 



- log(3y - 1) + i log(9y 2 + 3y + 1) 



l y 

5 5(-3^ + i) 



-3y 3 



dy= m 



— 4v / 3arctan ' 



1 - 27y 3 



V3 



1 

25' 
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where the constant — is found using (I25|) . 
25 

We plot the first two orders of the formal constant of motion in Fig. [3J 



Real part of the formal constant of motion 




-4.45 — 
-4.5 - 

-4.55 - 
-4.6 - 

-4.65 - 
-4.7 - 



60 80 100 120 140 



Imaginary part of the formal constant of motion 



20 40 60 80 100 120 140 

Figure 3. Formal constant of motion with Fq and F\. 

Since this formal constant of motion is almost a constant along any path in the 
same i?-domain, it can be used to find the solution asymptotically, writing 

5 x 
Placing the term log(3y — 1) (cf. (131))) ) in the equation above on the left side and 
C on the right side, taking the exponential, and solving for y, we obtain 



(39) y — i exp f—C—x+^logx + ( V3 — j arctan ^^L"*" ^ 



\ log( V + 3, + 1) + \ + k + )4 



The reason for taking the exponential in (|39[) is to take care of the branching 
due to logs, whereas the other log and arctan do not matter since the solution 
does not encircle their singularities. Equation (|39p contains, in an implicit form, 
the solution y to two orders in x. y can be determined from this implicit equation 
in a number of ways; we chose, for simplicity to numerically solve the implicit 
equation using Newton's method. The solution is plotted in Fig. @J where we 
take C = 2.18 — 4.65i and calculate the solution for the second half of the path 
corresponding to \x\ > 61.4. Note that the relative error is within 1.5%. 

Since the accuracy of the formal constant of motion is unaffected by going along 
the solution path as long as |a;| is large, we can obtain quantitative behavior of the 
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Real part of the solution 
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1 
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Constant of Motion 



0.2 
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80 
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0.5 
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Imaginary part of the solution 
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0.015 
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140 
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0.01 



0.005 
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Relative error 
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u 



60 80 100 120 140 



Figure 4. Comparison of solutions obtained numerically by the 
Runge-Kutta method and using the formal constant of motion . 



solution for very large |a;|. By contrast, in a numerical approach, the further one 
integrates along the path, the less accurate the calculated solution becomes. 



4.2. Finding the positions of the singularities. We illustrate how to find sin- 
gularities of the Abel's equation using Proposition UTI It is known [IT] that there 
are only square root singularities, and they appear in two arrays. 

For simplicity we choose a simple singular path along which y goes to +00. 

According to Proposition [Til we have 



(4(1, fi,(//) = / _ 3s l +1 J S 



arctan (^±i) + log(3y - 1) - \ log(9y 2 + 3y + 1) + ^ 



5 7 (-3y 3 + i) 2 

I .-27;, ; 1 _x tjau.i>uaJl(-^^ 



^— J- 2 v / 3arctan(22+l) + 21og(3y - 1) - log(9y 2 + 3y + 1) 
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1.3 



1.2 
1.1 



0.9 
0.8 
0.7 
0.6 



"i06.4 106.6 106.8 107 107.2 107.4 107.6 107.8 108 

FIGURE 5. A small section of the left-top plot in Fig. [4] 

Thus the position of the singularity is given by the formula 
(41) Xl =(7 + o(l) = x 

+ + ( l0g(3y ° 

where the initial condition (xo,yo) satisfies |cco| is large and yo is not close to any 
of the three roots. We note that the presence of the arctan in the leading order 
implies that the solutions remain quasi-periodic beyond the domain accessible to 
the methods in [TT]. In (|4T]) we have the freedom of choosing branch of log and 
arctan, which enables us to find arrays of singularities. 

For example, the position of a singularity corresponding to the initial condition 
.T = 10 + 60i, y = 0.7 + 0M, calculated using gU is x x = 9.80628 + 60.2167z, 
which is accurate with six significant digits, as checked numerically. 

The detailed behavior of the solution near the singularity can be found by ex- 
panding the right hand side of (|41[) . We omit the calculation here since there 
are many other methods to determine this behavior (cf. 17]) and it is of lesser 
importance to the paper. 

4.3. Connecting regions of transseries. We choose a path consisting of line 
segments The path in x consists of line segments connecting 50i, 50, — 50i, —50, 
50i, 50, — 50i, and — 50(-\/3 + i). This corresponds to an angle of 2n in the original 
variable, with initial condition y(50i) = 0.6. 

Along this path, the solution of ([35]) approaches all three complex cube roots of 
1/27. For instance, the root 1/3 is approached when x traverses the first quadrant 
along the first segment, the root (— 1) 4 / 3 /3 is approached when x goes to the lower 
half plane, and the root (— 1) 2 / 3 /3 is approached when x goes back to the upper 
half plane. Some of these values are approached more than once along the entire 



Runge-Kutta 




V3(l-^Harctan^ + 1 ^ * 
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path. This behavior can easily be shown using the phase portrait of Go, cf. (fT6|) 
and Corollary HI 

Note that along a straight line x — xq + xe t% where the angle t is fixed the leading 
term (with only Go on the right hand side) of the ODE (|35|) can be written as 

9, 



5(:eo + xe tl ) 
Denoting y± = Key and j/2 = Imy, we have 

= —3yf cos t — 3y| sin t + 9yiy| cos i + 9y??/2 sin £ 

ax 

-r^ = — 3y? sin £ + 3y| cos i + 9j/iy| sin i — 9y?y2 cos t 
ax 



cost 

~9~ 

sm t 

9" 



We can then analyze the phase portraits. For the purpose of illustration, we 
show some of them in Fig [5] and [7J 




-1.0- I \ \ 




FIGURE 6. Phase portrait of Re (y) and Im (y) for t = —tt/4. The 
" x " marks are the three roots. 



On the line segment connecting 50i and 50, it is clear that the initial condition 
0.6 is in the basin of attraction of 1/3 (cf. Fig. [5]). 

(-1) 4/3 

Since the only stable equilibrium is ao — 7; j on the line segment connecting 

o 

50 and — 50z the solution converges to ao (cf. Fig. [7]). 
Numerical calculations confirm this, (cf. Fig [5]). 

Note 14 (Absence of limit cycles). Finally, note that there cannot be a limit cycle in 
the phase portraits drawn if x goes along a straight line. If the solution y approaches 
a limit cycle, it must lie in an R-domain. Thus the formal constant motion formula 
(0) is valid, and the first term Fq specifies a direction for x. If x goes strictly along 
this direction towards 00 then the term alogx, which does not vanish in our case, 
will go to 00, contradicting the results about the constant of motion. On the other 
hand, if x goes in a different direction, then —x + Fo(y) goes to 00 much faster than 
alogx, again a contradiction. 
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, / , A \ \ 

/ / H \ V r X/ ^ j 



FIGURE 7. Phase portrait of Re (y) and Im (y) for t — 5ir/4. 




Length of the path from the starting point 

Figure 8. Behavior of the solution across transseries regions. 
Dotted horizontal lines are the imaginary parts of the three roots. 
The horizontal line is arclength. The path in x consists of line 
segments connecting 50z, 50, — 50z, —50, 50i, 50, — 50«, and 
— 50(v / 3 + i). This corresponds to an angle of 2tt in the original 
variable. 



4.4. Extension to higher orders. For higher orders, such as the Painleve equa- 
tions PI and P2, a similar procedure works, though the details are quite a bit more 
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complicated, and we leave them for a subsequent work. We illustrate, without 
proofs, the results for PI, y" = 6y 2 + z. Now, there are two asympotic con- 
stants of motion, as expected. The normal form we work with is u" + u'x~ 1 — 
u — u 2 /2 — 392x~ 4 /625 = 0. Denoting by s the "energy of elliptic functions" 
s = u' /2 — ti 3 /3 + u 2 (it turns out that s is one of the bicharacteristic variables 
of the sequence of now PDEs governing the terms of the expansion; thus the pair 
(it, s) is preferable to (u, u')), one constant of motion has the asymptotic form 

C\ = x — L(s, u) + x~ 1 Ki(s, u) + ■ ■ ■ 

In the above, denoting R = \/it 3 /3 + u 2 + s, L is an incomplete elliptic integral, 
L = j R^ 1 (s, u)du and the integration is following a path winding around the zeros 
of R. The functions K\, K2, ■ ■ ■ have similar but longer expressions. We note the 
absence of a term of the form a log x (the reason for this is easy to see once the 
calculation is performed). A second constant can now be obtained by reduction 
of order and applying the first order techniques, or better, by the "action-angle" 
approach described in the introduction. It is of the form 

C 2 = xJ(s) + [L(s)J(u, s) - J{s)L(u, s)] + x~ 1 K 1 H 

where J(u, s) — J R(s, u)du; when the variable u is missing from J(u, s) or R(u, s), 
this simply means that we are dealing with complete elliptic integrals. There is 
directionality in the asymptotics, as the loops encircling the singularities need to 
be rigidly chosen according to the asymptotic direction studied. A slightly different 
representation allows us to calculate the constants to all orders. Because of direc- 
tionality, a different asymptotic formula exists and is more useful for the "lateral 
connection" , that is, for calculating the solution along a circle of fixed but large 
radius, which will be detailed in a separate paper, as part of the Painleve project, 
see e.g. [9]. 
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